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, Abstract. We present a proposal for a gravitational wave detector, based on the 

' excitation of an electromagnetic mode in a resonance cavity. The mode is excited 

due to the interaction between a large amplitude electromagnetic mode and a quasi- 
I monochromatic gravitational wave. The minimum metric perturbation needed for 

' detection is estimated to the order 7 x 10^^'^ using current data on superconducting 

lO ■ niobium cavities. Using this value together with different standard models predicting 

\l ' the occurrence of merging neutron star or black hole binaries, the corresponding 

■ detection rate is estimated to 1-20 events per year, with a 'table top' cavity of a 

P\J , few meters length. 

o 
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During the last decades the quest for detecting gravitational waves has intensified. 
The efforts have been inspired by the indirect evidence for gravitational radiation 
[T], advances in technology and the prospects of obtaining new useful astrophysical 
information through the development of gravitational wave astronomy 25 . A number of 
ambitious detector projects are already in operation or being built all over the world, 
for example Ligo and ALLEGRO in USA, Virgo, AURIGA and GEO 600 in Europe, 
Tama 300 in Japan, AIGO and NIOBE in Australia 3J. Furthermore, there are well 
developed plans to use space based gravitational wave detectors, i.e., the LISA project 
|3]. The detection mechanisms are basically of a mechanical nature in the cases above, 
but there have also been several proposals for electromagnetic detection mechanisms [Hj. 

In the present paper we will investigate a detection mechanism based on the 
interaction of electromagnetic modes and gravitational radiation in a cavity with highly 
conducting walls. The main feature of our proposed gravitational wave detector 
is that it supports two electromagnetic eigenmodes with nearby eigenfrequencies, a 
possibility that has previously been discussed in Refs. jHl Ej • If one eigenmode is excited 
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initially (called the pump-mode), and a quasi- monochromatic gravitational wave with 
a frequency equal to the eigenmode frequency difference reaches our system, a new 
electromagnetic eigenmode can be excited due to the gravitational-electromagnetic wave 
interaction. The coupling mechanism is similar in principle to the wave interaction 
processes described in, e.g., Ref. jH]. 

The low-frequency nature of the gravitational modes, as compared to the 
electromagnetic resonance frequencies, at first seem to greatly limit the efficiency of 
a cavity with nearby frequencies. To get a large gravitationally induced mode-coupling 
for a simple cavity geometry, the estimated cavity dimensions become prohibitively 
large, i.e., comparable to the wavelength of the gravitational wave. Solutions to this 
problem has been found by Refs. [7j, who has considered a gravitational wave detector 
consisting of two coupled cavities. Cavities based on these principles have been built, 
and experimental results are presented in Refs. jH]- In this letter we consider a single 
cavity with a variable crossection. The main purpose of varying the crossection is 
the following. For a cavity with dimensions much smaller than the wavelength of the 
gravitational wave, the suggested geometry greatly magnifies the gravitational wave 
induced mode-coupling. In our present work we have simulated the effect of a varying 
crossection, by considering a cavity filled with three different dielectrics, in order to be 
able to perform most of the calculations analytically. It is straightforward to make a 
semi-quantitative translation of our results to the case of a vacuum cavity with a varying 
crossection. Using current data on the latter type of cavity [10111111121, we estimate the 
minimum detection level of the metric perturbation to the value /imin ~ 7 x 10~^^, where 
we have considered an inspiraling neutron star or black hole pair as a gravitational wave 
source. If such a level of sensitivity can be reached, neutron star or black hole binaries 
close to collapse could be detected at a distance r^ax ~ 10® ly. Adopting data from 
Ref. for the occurrence of compact binary mergers, we obtain the estimate 1-20 
detection events per year. 

In vacuum, a linearized gravitational wave can be represented by g^^ = rj^u + h^^iO 
where = x — ct, and h^x = —hyy = h^, h^y = hyx = in standard notation. 

Neglecting terms proportional to derivatives of and h^, the wave equation for 
the magnetic field is ^1 



and similarly for the electric field. Here n is the index of refraction. For the moment, 
we will neglect mechanical effects, i.e., effects which are associated with the varying 
coordinates of the walls due to the restoring forces of the cavity. | 

The coupling of two electromagnetic modes and a gravitational wave in a cavity will 
depend strongly on the geometry of the electromagnetic eigenfunctions. We can greatly 
magnify the coupling, as compared to a rectangular prism geometry, by varying the 

I This is true neglecting thermal noise and assuming a long cavity, as compared to the speed of sound 
over the gravitational wave frequency. 




(1) 
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cross-section of the cavity, or by filling the cavity partially with a dielectric medium. 
The former case is of more interest from a practical point of view, since a vacuum 
cavity implies better detector performance, but we will consider the latter case since it 
can be handled analytically. However, we will show how to make a semi-quantitative 
translation of our results to the case of a varying cavity cross-section. 

Specifically, we choose a rectangular cross-section (side lengths Xo and yo), and we 
divide the length of the cavity into three regions. Region 1 has length li (occupying 
the region < z < h) and a refractive index rii. Region 2 has length I2 (occupying the 
region li < z < li + 12), with a refractive index 712, while region 3 consists of vacuum 
and has length (occupying the region h + I2 < z < h + I2 + h = 0- The cavity 
is supposed to have positive coordinates, with one of the corners coinciding with the 
origin. Furthermore, we require that I2 <^ h, and that the wave number in region 2 
is less than in region 1. The reason for this arrangement is twofold. Firstly, we want 
to obtain a large coupling between the wave modes, and secondly we want an efficient 
filtering of the eigenmode with the lower frequency in region three. 

The first step is to analyze the linear eigenmodes in this system. The simplest 
modes are of the type 

iujxo ~ . ( m7rx\ 



Ey — Bzj sin 

mn \ xq 



j sm[kjZ + (Pj]e~"^\ (2a) 

TniTT \ 

B, = B,j cos ( sm[kjz + iPj]e-"^\ {2b) 



Xo J 

B^ = -^-^B^j sin ( ) cos[kjZ + ^j]e-'''\ (2c) 

in regions j = 1, 2 and 3, respectively, where the wave in region 3 is a standing wave, 
and m is the mode number. In region 3 we may also have a decaying wave 

Ey = -B,s sin sinh A;3^ + <^3 e"""*, 3a 

mTT V ^0 / 

B, = B^3 COS sinh[A;3-z + (^aje""^*, (36) 

\ Xo J 

rj hxo^ . f m7ix\ 

B^ = -B^asm cosh[A;3^ + ip-i\e (3c) 

V ^0 / 



mn 



Using standard boundary conditions, the wave numbers are calculated for an eigenmode, 
and the relation between the amplitudes in the three regions is found, and thereby the 
mode profile. We are interested in the shift from decaying to oscillatory behavior in 
region 3. Denote the highest frequency mode which is decaying in region 3 with index a, 
and the wave number and decay coefficient with kia, and k-^a respectively. Similarly, 
the next frequency, which is oscillatory in both regions, is denoted by index b. If we have 
/ ^ Xq (and m the same) these two frequencies will be very close, and a gravitational 
wave which has a frequency equal to the difference between the electromagnetic modes 
causes a small coupling between these modes. An example of two such eigenmodes is 
shown in fig. 1. 
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Figure 1. Mode-profiles for the eigenmodes are shown for the parameter values 
rii = 1.22, 71,2 — 1.0005, Iitti/xq = 143, l2'm/xo = 12.3 and l^m/xo = 184. i?norma is 
the solid line and Bnormb is the dotted line, and Ic — 0.27. 



We define the eigenmodes to have the form Bza,b = Babbit) Bnorma,b{z, y, z), where 
Ba^b is a time-dependent amplitude and the normalized eigenmodes -Bnorma.fe satisfy 
Jy \B^oTma,bf dV = V. We let all electromagnetic field components be of the form A = 
Aa{r) exp{—iujat) + Ah{r) exp{—iujht) + c.c. , where c.c. stands for complex conjugate, and 
the indices stand for the eigenmodes discussed above. The gravitational perturbation 
can be approximated by /i+_x = /i+,x exp{—iugt) + c.c, where we neglect the spatial 
dependency, since the gravitational wavelength is assumed to be much longer than all 
of the cavity dimensions. During a certain interval in time, the frequency matching 
condition Ub = uja + ojg will be approximately fulfilled. Given the wave equation ((l}, 
and the above ansatz we find after integrating over the length of the cavity 

^-^[ll-^)Bb = -hABah, (4) 

where 

^<^=Y^j^ ^^^5norm b dV, (5) 

and we have added a phenomenological linear damping term represented by 7. Thus 
we note that for the given geometry, only the /i+-polarization gives a mode-coupling. 
Calculations of the eigenmode parameters show that Ic may be different from zero when 
rii 2 7^ 1, and generally Ic of the order of unity can be obtained, see fig. 1 for an example. 
From Eq. (jlj), we find that the saturated value of the gravitationally excited mode is 

Bb... = ^hM^i^ , (6) 

In fig. 1 it is shown that we can get an appreciable mode-coupling constant 
Ic for a cavity filled with materials with different dielectric constants, and it is of 
interest whether or not this can be achieved in a vacuum cavity. As seen by Eq. 
(0), the coupling is essentially determined by the wave numbers of the modes, given 
by A; = (m^7r^/a;g — n^uj"^ / c^Y^'^ . Thus by adjusting the width xq in a vacuum cavity. 
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we may get the same variations in the wave numbers as when varying the index of 
refraction n. The translation of our results to a vacuum cavity with a varying width 
is not completely accurate, however. When varying xq, the mode-dependence on x and 
z does not exactly factorize, in particular close to the change in width. Moreover, the 
contribution to the coupling Ic in each section becomes proportional to the corresponding 
volume, and thereby also to the cross-section. However, since most of the contribution 
to the integral in Eq. (jSJ comes from region 1, our results can still be approximately 
translated to the case of a vacuum cavity, by varying xq instead of n such as to get the 
same wavenumber as in our above example. 

We denote the minimum detection level of the excited mode with Bhmm and 
the maximum allowed field in the cavity with -Bamax- The magnetic field -Bbmin is 
related to the minimum number of photons Ntm needed for detection by -B^min ~ 
(NbmfJ'O^b^/xoyol)^^'^ . Furthermore, we have ~ {l/3Q){uj'^/kj^c^) ~ 2/3Q, 

where Q = 07^/2717 is the quality factor of the cavity, and thus we obtain, using Eq. 

1 /2 

hmin ~ ; — ■ (7) 

Before giving estimates of the parameter values, we will investigate certain other effects 
that may limit the detection efficiency. 

Massive binaries produce monochromatic radiation to a good approximation, but 
close to merging, the frequency will be increasing rather rapidly which means that the 
phase matching will be lost. The cavity is designed to detect the frequency Ug, and we 
have assumed that the signal varies as exp(— icugt), but in reality we have exp{—iilj{t)), 
where for simplicity we assume ip = 0, dtp/dt = Ug at t = 0. The coherence time tcoh 
is roughly defined by cOgtcoh — 2-ip(tcoh/2) = vr. Thus, tcoh ~ [2tt/ (dug/dt)]^/'^ provided 
'^g^coh ^ 1- Using Newtonian calculations of two masses M in circular orbits around the 
center of mass, complemented by the quadrupole formula for gravitational radiation, we 
find 

dug _ 388^1 fGMugV^'' 

~dr-^^[-^) ■ 

Close to binary merging, the time of coherent interaction will be shorter than the photon 
life time, and in that regime the growth of the excited mode is limited by decoherence 
rather than the damping due to a finite quality factor of the cavity. However, formula 
^ can still be applied if we simply replace Q by Ubtcoh- 

To be able to estimate the number of photons needed for detection, we must study 
various sources of noise. The simplest effect is direct thermal excitation of photons 
in mode b. Since each mode has a thermal energy level of order KT, the number of 
such photons is ~ KT/hwb- However, we will also have a contribution associated with 
thermal variations in the length of the cavity. A standard model for the variations in 
length 61 is [6^ 

d'^61 uji d6l 2c-7 7 27 /• \ /^\ 
+ g +uJidl = -lujgh+ ey.Y>[iujgt) + ath, (9) 
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where ui is the eigenfrequency of longitudinal oscillations (of the order of the acoustic 
velocity in the cavity divided by the length), Qmec is the mechanical quality factor 
associated with these oscillations, and ath is the stochastic acceleration due to the 
thermal motion, giving rise to a random walk in the oscillation amplitude. First we 
note that the amplitude of the length variations with the gravitational frequency is 
6lg = Ugh+/{Ug — ujf). We will consider the case when ujf ~ 2a;^, which implies that 
the amplitude of the gravitational length perturbation is essentially unaffected by the 
restoring force of the cavity, as is the number of gravitationally generated photons. 

The thermal fluctuation contribution to the right hand side of Eq. ((3)), via the 
coupling to the pump wave, becomes proportional to exp[i(a;b — uja)t\. Since 
5lth ~ exp(zcij;t), the longitudinal oscillations give rise to slightly off-resonant (i.e., 
driven) fields with a frequency difference 5uo = uoi — uog compared to mode h. However, 
due to the stochastic changes in amplitude and/or phase of the longitudinal oscillation 
(where the time-scale for significant changes is given by r ~ Qmecl^i)i a contribution to 
mode h of the order {5lth) \/oJi/ 5ujQmecdt ~ KT / muj'f -^f 1 / Qmec {27r/6uj) is also made 
during a single oscillation period dt = 27i/6u!. Here m is the mass taking part in the 
longitudinal oscillation. During a time of the order tcoh, this contribution adds up to 
approximately (a/^vt/ 6uj) a/ KT/ mojf Sujtcoh/Qmec by a random walk process. Using 
5uj ^ ujg, the condition for the gravitational contribution ~ hltcoh to be larger than that 
of the thermal fluctuations can be written 



Assume that we want to reach a sensitivity hmin ~ 7 x 10~^^. We let xq = 0.1m, 
Uq = 0.2m, / = 4m , and m = 4 which gives Uf, ~ Attc/xq = 4 x 10"'^°rad/s. 
Furthermore, we take Ug = lOOOrad/s and let ui = 1500rad/s. We assume T = IK 
together with m = 2500kg. From ()10|) . we then find that we need a mechanical 
quahty factor Qmec ~ 6 x 10^ to reach the desired sensitivity. § Moreover, assuming 
the necessary number of photons for detection to be N^m = 20 (well above the direct 
electromagnetic noise level of a few photons), Ic ~ 0.3 and -Bamax = 0.2 T [TTJ [T2], 
we need an electromagnetic quality factor Q = 5 x 10^ to reach the desired sensitivity 
^min ~ 7 X 10~^^. Note that the coherence time is slightly longer than the photon 
life-time, as needed. 

Following the example given above for calculating the coherence time, i.e., a binary 
consisting of two compact objects, each of one solar mass Mq, separated by a distance d, 
the amplitude of the metric perturbation at a distance r is given by /i ~ 2G'^M'^ / (?rd. 
For definiteness we choose d = 80 km corresponding to Ug = lOOOrad/s, and thus 
for /iiiiin of the order of 7 x 10^^^ we obtain the maximum observational distance 
'^max ~ 8 X 10^ ly. Using data from Ref. we deduce that the number of galaxies 
within the observational distance rmax ~ 10*^ ly is of the order = 10^. Combining 
these figures with the expected number of compact binary mergers per galaxy and 

§ Note that mechanical quahty factors can be several orders of magnitude better, see e.g., Ref. 




(10) 
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million years [R]; we obtain a detection rate in the interval 1 — 20 events per year (the 
uncertainty is due to different models used for the birth of compact binaries). 

Our proposal for a gravitational wave detector is to a large extent based on currently 
available technology ^Sl UHl ^1 and our requirements are moderate given the 
performance of some existing microwave cavities. For example, values of Q ~ 4 x 10^° 
has been reported in Ref. JUl and quantum non-demolition measurements of single 
microwave photons have been made in Ref. 18j. Furthermore, the key performance 
parameters of the superconducting niobium cavities, i.e., the quality factor and the 
maximum allowed field strength before field emission, have been improving over the 
years, suggesting that the detection sensitivity can be increased even further. 

A sensitivity /imin ^ 10~^^ seems extremely good, but, on the other hand, 
idealizations have been made when making the estimate. In addition to any effects 
induced by the gravitational wave, the walls generally vibrate slightly due to the 
electromagnetic forces exerted by the large pump field. While these later oscillations 
clearly will be larger than the variations in length that are directly due to the 
gravitational wave, the associated nonlinearities will be harmonics of the pump 
frequency, and thus such effects do not couple to the other eigenmodes of the cavity. 
Furthermore, we have assumed that the detection of the excited mode is not much 
affected by the presence of the pump signals. Even though mode a is partially filtered 
out in region 3, the small frequency shift between the two electromagnetic modes may 
pose a certain difficulty in this respect: In particular, a very narrow bandwidth (pump) 
antenna signal must be used, in order to exclude the slightest initial perturbation at the 
gravitationally excited frequency Uj,. 

However, although there are technical difficulties in constructing an electromagnetic 
detector, the real advantage is the possibility to reduce the size of the devise. In the 
example presented above, the length of the cavity has been taken to be 4 m, and the 
cross section roughly 1/50 m^. This alone could prove to be useful when trying to set 
up new gravitational wave observatories. 
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